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A.  Description  of  the  Scientific  Research  Goals 

The  first  objective  of  this  research  is  to  develop  and  validate  a  novel  approach  to 


subgrid  scale  modeling  based  on  vector  wavelets.  This  approach  must  be  capable  of 
treating  non-equilibrium  turbulence.  The  second  objective  is  to  develop  a  LES  scheme, 
which  uses  wavelet-type  helical  basis  functions  for  the  space  resolved  scales,  to  describe 


anisotropic  and  wall-bounded  flows.  The  third  objective  is  to  integrate  the  wavelet  SSM 


in  this  new  LES  scheme.  After  developing  this  wavelet  LES  the  next  objective  is  to  carry 


out  LES  of  some  nonequilibrium  flows  which  will  lead  to  a  general  purpose  LES  scheme 


capable  of  providing  improved  estimates  of  dynamical  characteristics  such  as  skin 


friction  and  mixing,  as  well  as  enabling  prediction  and  adaptive  control  of  separation  in 


B.  Significant  Results  in  the  Past  Year 
Bl.  Introduction 

We  have  laid  the  foundation  to  develop  and  validate  LES  using  wavelets  as  a 
functional  basis  and  based  on  subgrid  scale  (SGS)  modeling  using  vector  wavelets. 
Wavelet  LES  (WLES)  consists  of  subgrid-scale  model  equations  (SSM)  and  space- 
resolved  model  (SRM)  equations  .  Using  a  vector  wavelet  decomposition  of  the  velocity 
field  a  simple  model  for  locally  isotropic  turbulence  has  been  derived  from  the  Navier- 
Stokes  equation.  This  model,  which  involves  no  empirical  or  adhoc  parameter, 
incorporates  nonlocal  inter-scale  interactions,  reveals  backscatter  and  can  be  applied  to 
represent  small-scale  turbulence  in  LES  schemes.  Stationary  solutions  of  the  model 
equations  produce  the  Kolmogorov  k'5/3  inertial  spectrum  and  the  k^  infra-red 
spectrum.  We  have  completed  derivation  of  the  SRM  equations  based  on  the  helical 
wave  decomposition.  We  will  test  the  SRM  equations  using  the  computational 
resources  in  this  NAS  operational  year.  A  wavelet-based  subgrid-scale  model  (WSSM) 
will  be  generalized  to  account  for  anisotropic  and  inhomogeneous  turbulence  in  wall- 
bounded  flows  and  will  employ  a  nonuniform  grid  to  resolve  the  near-wall  structures. 

B2.  LES  Scheme 

In  WLES  we  start  by  representing  the  flow  field  using  regularly 
distributed  fixed  helical  basis  functions  in  physical  space.  These  fixed  basis  functions 
correspond  to  the  macroscopic  (resolved)  scale  in  conventional  LES.  Then,  in  each  grid 
box,  we  consider  a  hierarchy  of  small-scale  wavelets  corresponding  to  the  subgrid  scales 
(SGS's).  Our  SSM  equations  are  derived  directly  from  the  incompressible  Navier-Stokes 
equations  by  a  Galerkin  projection  onto  a  complete  wavelet  basis.  Through  analytical 
approximations  (thereby  eliminating  calibration  constants)  to  the  wavelet  interaction 


2 


coefficients,  a  system  of  ODE's  is  obtained.  These  ODE's  describe  the  energy  contained 
in  wavelets  of  scales  (,  C/l,  /74,...,  //2E-1,  where  L  is  the  number  of  levels  in  the  SGS 
wavelet  hierarchy.  These  L  variables  enable  the  energy  spectrum  to  evolve  dynamically 
and  thus  make  the  SSM  ideally  suited  for  nonequilibrium  flows.  We  have  reproduced 
crucial  features  of  homogeneous  turbulence  using  this  model:  e.g.  infra-red 

spectrum,  the  value  of  Kolmogorov  constant  (1.48),  conservation  of  equilibrium  k^ 
spectrum,  self-similar  decay  of  stationary  spectrum  (with  the  correct  decay  rate)  etc. 
The  results  from  our  model  are  compared  with  recent  experimental  and  numerical^ 
results.  Currently  we  are  testing  the  model  for  nonequilibrium  situations.  We  are  also 
deriving  a  new  SSM  based  on  the  helical  wave  decomposition. 

For  the  WEES  scheme,  we  use  macroscopic  scales  resolved  on  Eulerian  grid 
points  (N  in  number)  while  the  residual  scales  are  modeled  using  wavelets  which  act  as 
a  crucial  dynamical  link  between  macroscopic  and  residual  scales.  This  feature  is 
particularly  suited  for  nonequilibrium  flows,  since  SGS  wavelets  mimic  the  physical 
nature  of  turbulence  and  allow  accurate  modeling  of  significant  spatial  and  temporal 
variations  in  turbulence  intensity.  A  notable  feature  of  this  scheme  is  that  most  of  the 
variables  are  used  to  describe  the  SGS  dynamics  (namely,  L  SGS  variables  for  each  grid 
point);  this  is  in  contrast  to  conventional  SGS  models,  which  have  no  parameters  for 
describing  the  dynamics. 

In  brief,  the  scheme  can  be  represented  in  the  following  way;  we  start  by 
decomposing  the  velocity  field  in  terms  of  the  wavelet  basis: 


where  Va  is  the  wavelet  corresponding  to  the  index  a  and  a  =  (N,v,m)  is  a 
composite  index.  N  represents  the  scale,  v  is  the  helicity  index,  and  m  is  the  location  of 
the  center  of  the  wavelet.  Aa  represents  the  amplitude  of  the  wavelets  used  to 
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decompose  the  velocity  field.  Note  that  N  =  0  represents  the  resolved  scales,  the  level  at 

which  boundary  conditions  are  specified.  To  decompose  velocity  field  at  the  resolved 

scales,  we  employ  a  pair  of  helical  basis  functions.  The  evolution  equation  for  the 

wavelet  amplitude  can  be  written  as 
d-A.(Y 

^=La(Ap) 


The  nonlinear  operator  L  consists  of  the  inviscid  and  the  viscous  part.  The  complete 
evolution  equation  can  be  written  as: 


yp  A=:yyT  aa+yk  a. 

^  Nvnjrm  N)J.m  NvnMjrmLXl  M|im  LA,1  ^  Nvnjun  Nirm 


This  equation  will  be  explained  in  detail  in  §B5. 

The  accuracy  of  the  scheme  depends  on  the  truncation  procedure  adopted  to 
calculate  the  coefficients  T  and  K  which  correspond  to  the  inviscid  and  the  viscid  part 
respectively.  The  influence  of  a  given  wavelet  decays  rapidly  as  we  move  away  from  its 
center.  The  truncation  procedure  is  based  on  a  threshold  approximation  which  retains 
all  the  coefficients  lesser  than  a  specified  value.  For  the  subgrid  scales,  the  interaction 
coefficients  are  assumed  to  be  1  if  a  triad  of  wavelets  overlap  and  are  set  to  zero  if 
otherwise.  All  calculations  are  carried  out  in  the  coefficient  space.  WSSM  has  no 
arbitrary  constants  and  thus  does  not  rely  on  experiments  for  empirical  data. 
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B3. Differences  with  the  traditional  LES 

Currently  popular  LES  methods  do  not  track  the  evolution  of  small-scale 
turbulence  and  thus  assume  that  at  each  instant  the  subgrid  turbulence  is  in  statistical 
equilibrium  with  the  large-scale  motion.  This  implies  that  all  scales  within  the 
nonstationary  part  of  the  spectrum  must  be  modeled  as  large-scale.  If  a  subgrid  model 
would  involve  variables  and  equations  for  describing  small-scale  turbulence  evolution, 
the  assumption  of  statistical  equilibrium  between  large  and  small  scales  would  no 
longer  be  required,  and  modeling  of  only  the  energy  containing  (inhomogeneous) 
largescales  would  be  sufficient.  This  will  significantly  reduce  the  cut-off  spatial 
frequency  for  the  large  scales,  thereby  decreasing  the  operation  count  of  the  LES 
drastically.  WLES  employs  just  such  a  model  using  a  vector  wavelet  approach  to  small- 
scale  turbulence  model. 

Another  significant  feature  of  the  WLES  is  that  it  does  not  have  a  sudden  loss  of 
resolution.  In  the  traditional  LES  schemes  above  the  cut-off  wave  number  suddenly  we 
have  only  a  statistical  description  of  the  flow;  this  accentuates  the  predictability 
problem  of  the  LES  drastically.  In  WLES  below  the  cut-off  wave  number  we  solve  a 
simplified  ODE  version  of  the  NSE  thus  preserving  a  model  dynamical  description  of 
the  flow.  This  method  is  inherently  consistent  with  the  fact  that  there  is  no  spectral  gap 
between  large  scale  structures  necessary  for  the  existence  of  turbulence  and  the 
secondary  eddies  generated  from  them  .  In  this  scheme  it  is  very  easy  to  satisfy  the 
"stringent"^  requirement  that  the  low  frequency  part  of  the  spectra  not  be  affected  by 
the  procedure  used  to  remove  unwanted  degrees  of  freedom  from  the  dynamical 
description  of  the  flow.  Moreover  this  approach  is  amenable  to  construction  of  models 
for  production  process  as  the  wavelets  closely  parralel  the  actual  eddies. 
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B4  Importance  of  the  model  equations 

Derivation  and  testing  of  the  model  equations  for  WSSM  is  a  central  part  of  this 
work  because  it  is  these  model  equations  which  enable  the  WLES  to  track  the  dynamical 
evolution.  However,  it  is  of  utmost  importance  to  note  that  in  the  process  we  also  have 
a  simplified  ODE  version  of  the  NSE  which  can  accurately  describe  turbulence. 
Discovery  of  these  equations  may  have  far-ranging  implications  for  turbulence 
modeling  beyond  the  confines  of  the  narrow  context  of  this  particular  LES 
development. 

Normally  turbulence  is  treated  without  taking  helicity  explicitly  into  account; 
however,  recently  it  has  been  shown  that  description  of  turbulence  in  terms  of  two 
helical  basis  functions  is  conducive  to  a  better  physical  understanding.  In  the 
framework  of  WLES  the  use  of  helical  basis  functions  allows  us  to  conserve  two  integral 
constraints:  energy  and  helicity. 

B5.  Deriving  the  model  equations 

To  obtain  a  finite  dimensional  approximation  of  NSE,  we  use  the  method  of  weighted 
residuals  with  wavelet  based  trial  function. 


ii(Nt)=  X  \t 

^  Nvn  Nvn 


Nvn 


The  divergence-free  vector  wavelets  v  are  defined  as 
^  Nvn 


(1) 


K  -1/2 


Here 


P^  =  (77c/9)2^^; 


cos(s)-cos(2s)^ 


s  =  k2%-x^^). 


(2) 


(3) 


The  index  N  denotes  the  scale  of  the  wavelet.  Fourier  transforms  of  wavelets  are 
localized  inside  spherical  shells  in  the  wavenumber  space  {re  2^  <  k  <re  2^^M-  Each 
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wavelet  can  be  considered  as  an  axisymmetric  vortex  with  its  axis  along  one  of  the  unit 

A 

vectors  e  ^  (Fig.  1).  The  wavelet  centers  are  randomly  distributed  in  the  physical 
space  with  a  density  which  is  chosen  to  produce  a  uniform  approximation  of  the 

velocity  field  in  the  N-th  wavenumber  shell.  The  vorticity  magnitude  of  each  vortex  of 

—1/3 

the  scale  N  is  significant  only  in  a  spherical  region  of  radius  p  with  the  center  located 

_3 

at  Xj^^.  We  choose  Pj^  as  the  volume  of  the  wavelet  localization  region  in  the  physical 
space. 

The  Galerkin  projection  of  NSE  on  the  basis  (2)  has  a  form 


yp  A  =yyT  a  a+yk  a 

Nvn|J.m  Njim  NvnM)j.mLXl  Mpm  LA,1  Nvn|un  N)im 


where 


P  =  fd^xjv  V 

Nvimm  J  L  Nvn  N|im  J 


T  =  \  [dYIv  (v  V)v  +v  (v  V)v 

NvnM^mlAl  2  J  \  Nvn  Mjam  0.1  Nvn  01  Mpm  J 

and  K  =  f  d\  |  v  Av  | . 

NvnjLim  J  t  Nvn  NfimJ 


(4) 

(5) 

(6) 

(7) 


The  matrix  elements  T  describe  the  nonlinear  interactions  between 

NvnM|im01 

V  ,  v  ,  and  v  .  Consequently,  they  have  significant  values  only  if  the 
Nvn  M|im  01  ^  ^  ^ 

localization  regions  of  the  individual  wavelets  forming  this  triad  overlap  in  physical 

space.  These  elements  scale  with  the  wavelet  scale  as 


T 

(N+P)vn(M+P)|J  m(L+P)  A1 


25P/2 


T 

NvnMpmOl 


(8) 
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Conservation  of  kinetic  energy  (due  to  the  nonlinear  terms  in  equation  (4))  is 
guaranteed  by  the  following  property  of  the  matrix  T 


T  +  T  +  T  =0. 

NvnM|imLAl  MiJ.mLX,lNvn  lAlNvnMiim 


We  represent  the  wavelet  amplitudes  as 

A  (t)  =  A  (t)a  (t),  (10) 

Nvn  N  ^  Nvn 

where  the  new  variables  A  ,  later  called  the  shell  variables,  are  linked  to  the  kinetic 

N 

energy  Ej^(t)  in  the  N-th  wavenumber  shell  as 

E^(t)=  |pn  A^(t).  (11) 

Taking  an  ensemble  average  over  the  realizations  with  the  same  A^(t)  we  obtain  the 
following  exact  equations  for  A^^  from  (4),  (10),  and  (11) 


A  =  y<I>  AA+AA, 
N  ^  NML  ML  N  isf 


where 


<5  .  SIT 

NML 


vjLA  nml 


(a  a  a 

NvnM|imlAl  'Nvn  Mjim  L^l 


^  P>"I  I  K 

v)i  nm 


a  a 


~  ^^N  ^  ^  Nvnjim  '  Nvn  Nfxm 


Here  <>  denotes  the  statistical  ensemble  average  over  all  possible  realizations  of 

the  phase  variable:  a  (t). 

T  Nvn 
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To  obtain  the  model  equations  for  we  introduce  the  following  closure 


correlations 


^  Nvn  Miim  LXl 


(OT  +QT  +QT,  ),  (15) 

'  iSIML  NvnM|i,mlAl  MLN  M|imLA.]Nvn  LNM  LX]NvnM|im 


where 


^MML 


=  2 


-5N/2  -1/2 

^NML 


(16) 


and  Z.  „  „  is  the  number  of  the  overlapping  wavelet  triads  that  include  the  wavelet 

V  as  one  of  the  interacting  wavelets.  Note  that  /  a  a  a  \  is  scale  invariant 

Nvn  '  Nvn  M|xm  LAl  ' 

-1/2 

due  to  the  first  multiplier  in  (16),  while  Zj.^^  accounts  for  the  resulting  interaction 
between  the  wavelet  v  and  Z.  _  „  pairs  of  wavelet,  which  have  localization  regions 

Nvn  rsiJViL  ^ 

overlapping  with  the  localization  region  of  In  addition,  we  neglect  the  cross¬ 

correlations  of  the  normalized  wavelet  amplitudes  inside  the  shells 

/a  a  (17) 

The  factor  1/3  follows  from  (10)  and  (12).  This  assumption  is  necessary  only  to  simplify 
calculations  of  K  and  is  justified  because  the  value  of  nondiagonal  elements  are  much 
less  than  the  diagonal  elements. 

For  a  random  distribution  of  the  vortex  centers  the  sums  in  (12)  were  calculated 
analytically.  The  signs  of  are  uniquely  determined  by  energy  conservation  in  the 

limit  of  zero  viscosity.  The  diagonal  interaction  coefficients  (f.e.  when  two  of  the  three 

indices  are  equal)  <E)  ,0  ,  and  <I>  have  significantly  higher  values  than  the 

^  NMM  NMN  NNL 
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adjacent  nondiagonal  elements.  To  simplify  the  model  equations  and  make  them 
amenable  to  analytical  investigation  we  first  neglect  the  nondiagonal  coefficients  O 

NML 

and  then  approximate  the  diagonal  coefficients  by  exponential  functions.  This  leads  to 
the  following  system  of  model  equations 


A  =  2' 
N 


5N/2 


s 

M 


B 


M 


2'^  A  A 

N  N-M 


<3M/2  2 

^  ^+M 


2N 

-  2  DA 


N 


(18) 


where  B  =  0;  B=  0.387;  B=  0;  B  =  2.19;  D  =  937t^/35  .  Note  that  all 
M<-1  -1  0  M>1 


these  parameters  are  analytically  calculated  and  not  empirically  determined.  We 
choose  the  length,  velocity  and  time  scales  to  be  1^  =  ,  1^  /v  and  v/  1^  respectively, 

where  v  is  the  kinematic  viscosity.  In  these  units  it  is  convenient  to  consider  Reynolds 

number  for  each  scale  Re  =  A  2^^ 

N  0 


At  zero  viscosity  (the  same  as  putting  D  =  0)  the  equations  (18)  conserve  the  total 
kinetic  energy  E=^  E^and  has  the  exact  stationary  solution 


A 


N 


:  0.557e^^^ 


(19) 


9 /Q  /o 

which  corresponds  to  the  Kolmogorov  energy  spectrum  E(k)  =  1.478'  k^  ^  .  If  the 
energy  is  initially  confined  in  only  one  shell,  e.g.  at  t=0,  =  0  orily  if  N  0,  then  small 

•  N  2  . 

disturbances  SA^  grow  as  SA^  =  2.19-2  A^.  This  analytical  result  corresponds  to  the 


elevation  of  the  k^  infra-red  spectrum,  which  is  simultaneously  generated  at  all  scales. 

It  is  to  be  noted  that  the  closure  relation  (15)  was  derived  later.  Initial  derivation 
of  the  equation  (18)  followed  dynamical  assumptions  to  approximate  the  interaction 
coefficients.  The  approximation  procedure  for  these  interaction  coefficients  is  presented 
in  the  appendix  1. 
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B6  Results  from  the  model  set  of  equations 

Solutions  of  (18)  were  obtained  numerically  for  viscous  cases  by  a  variable 

coefficient  ODE  solver^  with  a  pre-conditioned  Krylov  subspace  iterative  method^.  The 

most  significant  result  so  far  is  presented  in  the  Fig  2.  The  initial  energy  was  confined 
in  the  shell  N  =  0  and  the  initial  (Re^_^)  was  lO^.  The  solutions  display  a  infrared 

spectrum  to  the  left  of  the  initial  energy  containing  scale  and  k‘^/3  Kolmogorov  inertial 

spectrum  to  the  right.  Thus  both  the  energy  transfer  from  the  small  scale  to  the  large 

scale  (backscatter)  and  the  large  scale  to  the  small  scale  (cascade)  can  be  adequately 

captured  by  the  simple  model.  The  evolution  is  shown  from  0.1  to  1  (in  steps  of  0.1) 

eddy  turnover  time.  In  the  dissipative  range  the  energy  spectrum  decays  as  exp(-P'qk) 

with  p  =  5.37,  which  is  close  to  P=5.2  reported  recently  by  Saddoughi  and  Veeravallil  for 
a  large  Re  boundary  layer(Reg  =  74000;  R^  =  600).  The  continuous  spectrum  was 

obtained  by  a  cubic-spline  interpolation  of  the  energy  spectrum  corresponding  to  the 
shell  variables.  Our  results  are  compared  with  the  time  evolution  of  the  kinetic-energy 
spectrum  in  an  unforced  pseudo-spectral  128^  large-eddy  simulation  by  Lesieur  and 
Rogallo^,  started  initially  with  a  k^  spectrum  at  low  k.  The  point  to  note  here  is  that  the 
wavelet  based  model  equations  capture  the  infra-red  spectrum  with  the  same  accuracy 
as  an  LES  where  the  eddy  viscosity  was  evaluated  using  EDQNM  kinetic  energy 
transfer.  However,  our  simulations  cover  a  much  wider  range  of  scales  (  6  decades 
compared  with  2  decades  in  the  LES)  and  also  include  the  dissipation  range.  A  detailed 
comparison  between  the  LES  results  and  ours  is  beyond  the  scope  of  this  letter.  Note 
that  this  evolution  is  tracked  by  using  only  32  shell  variables  compared  with  128^ 
variables  in  the  LES. 

Thus  we  have  developed  a  successful  model  for  the  small  scale  turbulence 
dynamics  involving  simple  ordinary  differential  equations  which  are  much  simpler 
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than  the  EDQNM  model  (only  one  known  to  represent  backscatter  accurately)  and  thus 
have  a  good  chance  of  being  successful!  as  a  subgrid  model.  More  importantly,  since 
these  equations  can  track  the  small-scale  evolution  using  very  few  shell  variables  and 
no  adhoc  parameters,  they  can  be  used  to  model  nonequlibrium  turbulence. 

The  growth  and  decay  of  the  shell  variables  are  shown  in  the  figure  3.  It  is  to  be 
noted  that  the  there  is  a  delay  time  before  the  amplitudes  start  growing  which  increases 
as  the  wavenumber  increases.  This  delay  represents  the  time  required  for  the  energy  to 
cascade  down  to  smaller  scales.  After  an  initial  rise  in  the  amplitudes  the  value  of  the 
shell  variables  settle  down  to  the  value  of  the  Kolmogorov  constant.  As  the 
wavenumbers  grow  the  effect  of  viscosity  becomes  more  and  more  pronounced  in 
lowering  the  equilibrium  value. 

The  value  of  Kolmogorov  constant  was  estimated  from  the  rate  of  energy 

estimate  into  the  system.  The  spectrum  is  given  by  E(k)  =  Ce^/^k-^/S,  while  it  can  be 
shown  that  the  energy  of  the  n-th  shell:  E^^  ~  4.721e^^^A^^.  From  these  two  relations 

Kolmogorov  constant  was  estimated  to  be  1.48. 

In  figure  4  we  observe  the  free-decay  of  turbulence  when  the  forcing  is 
withdrawn.  The  spectrum  retains  the  self-similarity  as  expected.  The  rate  of  energy 
decay  was  proportional  to  t'^. 

In  figure  5  we  observe  the  conservation  of  an  initial  k^  energy  spectrum  to  the  left 
of  the  energy-containing  wavenumber.  This  is  a  further  test  to  observe  if  the  behavior 
of  the  equations  closely  parallel  the  behavior  of  turbulence.  The  k^  energy  spectrum 
correspond  to  the  situation  when  the  energy  is  equally  divided  among  all  scales  and 
hence  this  is  a  situation  where  there  is  no  flow  of  energy  into  these  scales.  This  feature 
was  correctly  captured  by  the  model  equations. 
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In  figure  6  we  test  the  behavior  of  the  equations  under  a  simple  case  of 
nonequilibrium  turbulence  where  we  add  a  periodic  driving  force.  We  observe  that  the 
flux  to  the  large  scale  is  quite  different  from  the  flux  to  the  small  scales  while  the 
dissipation  is  essentially  constant.  This  demonstrates  the  capability  of  the  model  to 
handle  non-equilibrium  situation.  This  extension  to  non-axisymmetry  will  be  the 
primary  focus  of  our  efforts  next  year. 

B7  Extension  to  helical  models:  gateway  to  anisotropic  and  inhomogenous  flows 

To  derive  a  helical  version  of  the  WSSM  it  is  not  advisable  to  proceed  directly  by 
using  wavelet  basis  functions.  It  is  better  to  work  first  with  a  helical  Fourier 
representation  of  the  Navier  Stokes  equation  and  then  introduce  wavelets.  This  not 
only  makes  the  analytical  calculations  tractable  but  also  allows  us  to  generalize  the 
WSSM  to  inhomogeneous  helical  turbulence  with  random  anisotropy.  This 
development  of  an  explicit  helical  Fourier  representation  of  NSE  is  succinctly  sketched 
below. 

In  the  helical  Fourier  representation  the  equations  are  equivalent  to  the  original 

NSE,  however  they  describe  the  motion  of  the  fluid  in  terms  of  four  scalar  functions:  the 
amplitudes  AQ^(k,t)  and  phases  0^(k,t)  of  the  right  handed  (a  =  -(-  1)  and  the  left 

handed  (a  =  -l)  helical  Fourier  modes.  It  is  possible  to  express  the  spectral  densities  of 
kinetic  energy  E  and  helicity  using  AQ^(k,t)  alone. 

E(k,t)  =  2  X  A2(k,t);  H(k,t)  =X  akAlikA). 

Thus  both  fundamental  quantities  E  and  H  have  diagonal  form  in  this  representation. 

The  equations  for  the  amplitudes  and  phases  of  the  helical  Fourier  modes  follow. 

0,+  kV„(k,t)  = 
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=  E  E  jd^q  M^^(k,q,p)cos[©^(k,q,p)-<l>Jk,t)  -O^Cq,!)  -<I>^p,t)]A|5(q,t)  A^p,t). 
A„(k,t)3,  4>„(k,t)  = 

=  Eljd^q  M„p,(k,q,p)  sin[0^p^(k,q,p)  -Ojk,t)  -<I>p(q,t)  -0^p,t)]Ap(q,t)  A^p,t), 
where  M^p^(k,q,p)  and  0^p^(k,q,p)  represent  the  anaplitude  and  the  phase  of  the 

kernel  of  the  integral  operator  describing  the  interactions  of  the  triads  of  helical  Fourier 
modes. 

The  next  step  is  to  derive  a  closure  procedure  for  the  helical  equations 
independent  of  the  wavelet  formalism.  Then  it  will  be  very  easy  to  extend  this  scheme 
to  anisotropic  and  inhomogeneous  flows. 

B8.  Conclusions 

We  have  developed  a  successful  model  for  the  small  scale  turbulence  dynamics 
involving  simple  ordinary  differential  equations  which  are  much  simpler  than  the 
EDQNM  model  (only  one  known  to  represent  backscatter  accurately)  and  thus  have  a 
good  chance  of  being  successful!  as  a  subgrid  model.  More  importantly,  since  these 
equations  can  track  the  small-scale  evolution  using  very  few  shell  variables  and  no 
adhoc  parameters,  they  can  be  used  to  model  nonequlibrium  turbulence. 

The  preliminary  WSSM  based  on  energy  amplitudes  predict  the  energy  spectrum 
quite  accurately.  Thus,  it  is  worthwhile  to  extend  this  WSSM  based  only  on  energy 
amplitudes  to  non-equilibrium  situations,  while  continuing  to  develop  helical  WSSM 
which  uses  both  energy  and  helicity  amplitudes.  Since  the  WSSM  that  uses  only  energy 
amplitude  is  much  simpler  than  the  helical  SSM  the  use  of  it  will  be  justifiable  if  it  can 
deliver  comparable  performance  to  helical  WSSM  even  if  the  domain  of  its  applicability 
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turn  out  to  be  somewhat  restricted.  Moreover  the  development  of  only  energy  based 
WSSM  will  provide  important  guidelines  and  checkpoints  for  the  helical  WSSM. 

Direct  attempts  to  derive  helical  WSSM  require  extremely  lengthy  analytical 
calculations;  unfortunately  these  calculations  are  not  amenable  to  solutions  via  standard 
analytical  softwares  either.  Hence  the  process  was  divided  into  three  parts:  i)  derivation 
of  NSE  in  the  Fourier  space  using  energy  and  helicity  amplitudes  ii)  implementation  of 
a  closure  procedure  and  iii)  introduction  of  wavelet  basis.  Derivation  of  the  equations 
are  complete  and  rest  of  the  work  is  in  progress.  Since  the  closure  procedure  is 
introduced  prior  to  the  introduction  of  wavelet  generalization  to  anisotropic  and 
inhomogeneous  turbulence  becomes  simpler  to  achieve. 

Appendix 


Using  the  above  relations 


f  X-  9371^ 

K.,  ,,  =  V,,  .  (k)(27tik|)2v^^  .  (k)  =  =  26.2 

NvnMjdm  J  Nvm  ^  ^  Mp.mi  ^  '  oD 


NvnM|imLX.L 


=  /  dxj  dk')J  dk")J  dk- 


-3 

3i  Y  y(N+M+L)  (k'x  ev)i  (k"  xev)j 


(k'"  X  ev)i  -  27ri  (k'*x  +  k"  •x  +  k"’  -x  -  k'»x  -  k"*x  +  k  -x) 

27tik'"^^i - ^e 

k"' 


The  exponential  part  can  be  rewritten  as 

-27ii(k'*x^,  +k"*x^^  +k"'*x  -(k'+k"+k"')*x) 

^  ^  Nn  Mm  LI 

Since  x  appears  only  in  the  exponential  we  integrate  over  x  first. 
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=  5(k'  +  ’k"  +  k"’). 


-27ti(k’+k"+k"')*x 

Now  e 

Hence  the  integral  exists  only  if  k'  +  k"  +  k'"  =  0  k'"  =  -k'  -  k" 

Changing  notation, 

k'  =  k  and  k'"  =  k'  ^  k"  =  -k’  -k . 

Thus  the  structure  of  the  R  matrix  is 

Rv.  ..  r.r  =  f  dk’ f  dk"fv^^  .(kW  .(-k'-k)2rak’V^.,.  (-k')|. 

NvnMumLXL  J  J  \  Nvni  ^  Mumj  ^  lAli  '  ') 


Now  we  proceed  to  obtain  the  model  equations  by  replacing  the  exact 
expressions  for  the  interaction  coefficients,  with  a  threshold  approximation.  The 
approximation  is  equivalent  to  stating  that  the  modelled  R  is  equal  to  the  rms 


value  of  exact  R  when  three  wavelets  overlap  and  0  when  they  do  not. 


triads 


Non  overlapping  wavelet  triad 
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For  some  fixed  levels  N,M,L  we  estimate  by  integrating  over  all  possible  configurations 


’<Nn  '  ’‘Mm  '  ’‘u- 


Rnv  nM  umL  A1 


This  involves  18  integrations:  3  for  dx^^^,  3  for  6  +  6  for  each  R  in  ,  and  27 
summations  over  v,  p  and  X. 


Actual  distribution  of  R  is  complex  so  it  is  approximated  by  a  step  function,  i.e.  if  all 
three  wavelets  (N,M,L)  overlap  in  some  region  R  has  some  constant  value;  else,  it  is 


zero. 


This  yields  W  = 


q  =  M+L  -N,  M  <N,  L<N 
q  =  M-N,  M<  N, 


Thus  a  step  approximation  of  R  as  described  above  yields  the  scaling  law  k‘(5/3). 

A  more  accurate  approximation  of  R  in  the  mean  square  sense  yields  +  X. 

X  is  a  small  correction.  In  the  above  estimate  of  R  the  sign  of  R  is  ambiguous. 

By  substituting  the  expressions  for  basis  functions  and  proceeding  as  before  it  can  be 
shown  that 

X  j  j  RNvnMHmLXldXj^^  dx^^j  =  2  ^(N+M+L)  J 

vjiA. 

Thus, 


k^k'^  (1  +  cos^9)sin^e' 
k^+k'^-2kk'cos0 


k^k'^  (1  +  cos^e)sin^9 
k^+k'^-2kk’cose  j 
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Point  to  note: 


We  can  write  =  5 

This  will  make  the  coefficient  matrix  symmetric. 

However,  in  this  case, 

1  2  1  f  2 

Tnml  =  4  Rnml  +  2  Alm  ^  4 

Thus  we  need  to  evaluate  the  cross  terms,  which  can  be  done  in  the  following  manner: 


VjiX 

1 


[vnMumLW 


32^  3(N+m+l)  f  (k-k'cos9)sin^9cos9 

73  J  1  k^+k'^-2kk'cose 


2'^-'<k<2^ 
2^'^  <  k  <  2^ 


2^’^  <  A/k^T^-^kk^o^  <  2^ 

Based  on  approximations  of  Tnml  we  can  arrive  at  the  model  equation 
'^NM  M  ”  ^tT\j ^N’ 

Mt^O 

Using  scaling  arguments  we  get  the  equation  in  the  form 


2  t 

=  2.19  2^  [2^  (An-1  -  ^n^n+I^  +  X  2  "  ^^hAj-p 

P=1 


13 
6  P 


Here  p  denotes  the  levels  of  interaction. 


C.  Plans  for  Next  Year’s  Research 

Our  first  objective  is  to  further  validate  our  wavelet  LES  scheme  by  comparing 
turbulence  statistics  with  DNS  data  for  turbulent  boundary  layers  and  a  backward 
facing  step,  at  low  Rx  and  at  equilibrium.  We  are  particularly  interested  in  capturing 
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prominent  features  of  these  flows,  such  as  the  near-wall  streak  spacing  and  the 
boundary  layer  reattachment  length  for  the  backward  facing  step. 

To  evaluate  the  applicability  of  wavelet  LES  for  studying  flow  physics  in 
practical  flows,  we  will  investigate  nonequilibrium  effects  such  as  strong 
inhomogeneity,  pressure  gradients  and  wall  curvature.  Results  from  wavelet  LES  for 
moderate  will  first  be  compared  with  those  from  other  LES  subgrid-scale  models, 
followed  by  wavelet  LES  at  R^  unattainable  by  conventional  schemes.  These  high  Rx 
results  will  be  validated  and  compared  with  experiments  on  turbulent  boundary  layers 
to  be  carried  out  concurrently  in  our  laboratory.  Also,  to  investigate  the  performance  of 
wavelet  LES  in  bounded  flows,  we  will  simulate  transition  to  fully  developed  turbulent 
pipe  flow  and  an  oscillating  pipe  flow  and  compare  these  results  with  available 
experimental  data. 

Validation  of  Wavelet  LES.  The  wavelet  LES  will  be  validated  by  comparing  with  DNS 
of  turbulent  charmel  flow  at  Re  =  33006  and  DNS  of  flow  over  a  backward  facing  step 7 
at  Re  =  6000.  For  the  backward  facing  step,  we  will  run  two  cases  corresponding  to  the 
low  step  (W/H  =  2.5)  and  high  step  (W/H  =  1.25),  where  H  and  W  are  respectively  the 
step  and  channel  heights. 

Investigation  of  nonequilihriunt  effects.  Wavelet  LES  can  in  principle  be  accurately 
applied  to  high  Rx,  in  the  range  Rx  =  10^  to  Rx  =  10^.  These  high  Rx  values  will  be  used 
for  experimental  comparisons  (data  from  turbulent  boundary  layer  facility  in  our 
laboratory)  and  studies  of  fundamental  turbulence  physics;  low  to  moderate  Rx  are 
necessary  for  comparisons  with  DNS.  First,  we  will  consider  fully  developed  turbulent 
boundary  layer  with  periodic  boundary  conditions  and  an  impulsively  applied  mean 
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spanwise  pressure  gradient  to  investigate  inhomogeneity  and  high  Re  effects.  We  will 
then  simulate  a  boundary  layer  with  spanwise  wall  motion,  for  comparison  with  work 
in  progress  by  Kim  et  al  To  investigate  the  effects  of  wall  curvature,  we  will  consider 
a  turbulent  flow  over  a  cylinder  for  y  =  5,  where  y  is  the  ratio  of  the  boundary-layer 
thickness  to  cylinder  radius.  These  results  will  be  compared  with  the  minimal  flow 
DNS  of  Neves  et  al.  To  investigate  resolution  of  separation  and  reattachment,  we 
will  perform  LES  of  flow  over  a  backward  facing  step  for  W/H  =  1.25  (high  step)  and 
2.5  (low  step),  with  no-slip  conditions  on  the  solid  wall. 

LES  of  bounded  flows.  For  comparison  with  experimental  data,  we  will  first  carry  out 
wavelet  LES  of  transitioning  pipe  flow,  for  comparison  with  DNS  and  experiments  by 
Eggels  et  al.  10.  The  Reynolds  number  based  on  the  centerline  velocity  and  pipe 
diameter  (Rec)  will  be  7000  in  this  case.  To  evaluate  our  LES  scheme  under  non¬ 
equilibrium  conditions,  we  will  consider  an  oscillating  pipe  flow  with  Rec  ~  10000. 
Experimental  data  for  this  flow  will  be  obtained  from  the  existing  pipe  flow  facility  here 
at  the  University  of  Houston. 

Extension  to  anisotropy  and  inhomogeneous  flows 

The  first  step  in  this  direction  will  be  to  derive  a  closure  procedure  for  the  helical 
equations  and  compare  the  results  from  this  closed  helical  SSM  with  the  WSSM  results. 
As  the  helical  SSM  introduces  wavelets  after  the  closure  procedure  it  can  easily  be 
tailored  to  be  applied  to  anisotropic  and  inhomogeneous  flows. 
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Status  compared  to  original  milestgnes 


Originally  proposed:  YEAR  I 

Achieved:  YEAR  I 

Repeat  tests  of  analytical  and  compu¬ 
tational  analysis  of  prior  hierarchical 

SGS  models. 

A  detailed  analysis  was  done  and  a 

way  of  improving  earlier  models 

was  found  which  led  to  a  new  model. 

A  priori  testing  of  basic  SGS  model 

using  DNS  data  from  homogeneous 

turbulence  and  free  shear  flows. 

A  new  SGS  model  was  derived  and 

tested. 

Develop  and  validate  the  basic  LES 

scheme  for  homogeneous  turbulence: 

to  be  started  in  first  year  and 

completed  in  the  second  year. 

The  basic  equations  for  LES  has  been 

derived  and  tested.  The  code  will  be 

developed  and  validated  this  year. 

Originally  proposed:  YEAR  II  &  III 

Plans:  YEAR  II 

Validate  LES  for  wall-bounded  flows 

against  DNS 

To  be  done  (Details  above) 

Refine  new  SGS  to  handle  anisotropy 

and  SGS  helicity 

To  be  done  (Details  above) 

Develop  SGS  for  nonequilibrium 

turbulence 

To  be  done  (Details  above) 
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E.  List  of  Figures 

Figure  1.  A  vector  wavelet  with  z  axis  as  the  axis  of  rotation.  The  hatched  section  shows 
the  velocity  profile  of  this  wavelet  in  the  x  direction. 

Figure  2.  Kinetic  energy  spectrum  of  unforced  freely-decaying  locally  isotropic 
turbulence,  (a)  results  from  wavelet  based  model  equations  display  the  infra-red 
range,  k'^/^  Kolmogorov  inertial  range  and  the  dissipation  range.  The  region 
enclosed  by  the  dashed  lines  display  the  corresponding  range  of  scales  for  the  LES 
shown  in  (b)  which  reproduces  1283  large  eddy  simulation  of  Lesieur  and  Rogallo^. 

Figure  3.  Growth  and  decay  of  the  shell  variables. 

Figure  4.  Free  decay  of  turbulence. 

Figure  5.  Conservation  of  initial  k^  energy  spectrum. 

Figure  6  (a).  The  nature  of  excitation 

Figure  6(b).  Comparison  between  the  energy  flux  from  and  toward  large-scales. 
Amplitudes  of  periodic  force  and  dissipation  are  not  plotted  to  scale. 
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